Few-Photon Electron-Positron Pair Creation by Relativistic Muon Impact on Intense 

Laser Beams 
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Electron-positron pair production in combined laser and Coulomb fields is studied. To this end, 
the Feynman diagram for multiphoton pair creation by muon impact on a circularly polarized high- 
frequency laser beam is evaluated within the framework of laser-dressed quantum electrodynamics 
employing relativistic Volkov states. In the limit of low laser intensity, the result is shown to coincide 
with the known expression for multiphoton pair creation by a proton which is treated as an external 
Coulomb field. A scaling of the total pair creation rate is analyzed. The recoil distribution is 
calculated numerically and its dependence on the projectile mass is discussed. Energy spectra of 
the created particles and angular spectra of the scattered muon are presented. 

PACS numbers: 12.20.Ds, 13.40.-f, 32.80.Wr, 42.55.Vc 



I. INTRODUCTION 



The creation of matter-antimatter particle pairs in ex- 
ternal electromagnetic fields is a characteristic effect in 
relativistic quantum theory. For example, a very strong 
static electric field may cause the spontaneous creation 
of e + e~ pairs, if the field strength surmounts the critical 
value of E cr = m^c 3 /eh — 1.3 x 10 16 V/cm found already 
by Sauter 1]. Here, m e denotes the electron mass, e 
the elementary charge unit, c the light velocity in vac- 
uum, and h Planck's constant. In recent years, interest 
in nonlinear pair creation processes via multiphoton ab- 
sorption from external laser fields has raised 0, B 0] ■ Yet 
a periodically changing electromagnetic plane-wave field, 
like a laser beam, cannot by itself lead to pair creation, 
independently of its frequency or intensity as shown by 
Schwinger [3]. Thus, an additional source is needed 
such as e. g. a second laser beam [f| 0, H, H E] , a single 
(non-laser) photon [HHGSEi, 

or a charged particle. 
Here, we consider e + e~ pair production by the collision 
of charged particles with laser fields. 
There are two possible channels via which the creation of 
the pairs can take place. The first channel is of Bcthc- 
Heitler type where the pair is created by a virtual pho- 
ton from the Coulomb field of the projectile particle and 
r real photons from the laser field. While Bethe and 
Heitler originally treated the linear case of pair pro- 
duction by a single high-energy photon (r = 1) [l5j], be- 
low we focus on the nonlinear Bethe Heitler effect which 
involves the absorption of r > 1 laser photons and thus 
depends on the laser intensity in a nonlinear way. In 
the case of a muon projectile, the symbolic equation for 
this process is fi + raj — > fi + e + e~ , with the laser fre- 
quency ui. In the second possible channel, the pair cre- 
ation takes place by the collision of the laser beam with a 
real photon (nonlinear Breit- Wheeler process [l||), with 
the latter stemming from a Compton scattering event, 
for instance. This indirect mechanism is important for 
light projectile particles, such as electrons. 
In 1997, laser-induced e + e~ pair creation has been ex- 
perimentally observed at SLAC (Stanford, USA). In this 



pioneering experiment, highly relativistic electrons with 
an initial energy of 46GeV were scattered by an optical 
laser pulse with an intensity of 10 18 W/cm 2 (l7l. H^|. The 
electrons were Compton back-scattered, thus emitting 
30GeV 7-photons which in turn collided with the laser 
beam. In this collision, r = 5 laser photons of Hui w 2eV 
combined their energies with the 7-photon and thus pro- 
duced an e + e~ pair via the nonlinear Breit-Wheeler pro- 
cess uj(j -\- tuj — ► e+e" [m, E0SG3. 
Inspired by the SLAC experiment, several theoretical 
studies bear on the — experimentally not yet observed — 
Bethe-Heitler creation of electron-positron pairs in the 
collision of a relativistic pro ton or nucleus with a laser 

beam [H E, E E Ell E E E E E H3- In 

the rest frame of the projectile particle, the energy of 
the laser photons is significantly higher than in the lab- 
oratory frame because of the relativistic Dopplcr shift. 
Recently, e + e _ pair creation via the collision of a high- 
energy neutrino with an intense laser beam was calcu- 
lated by TlNSLEY [3l[. There, the creation of the e + e~ 
pair takes place via an intermediate Z° boson. 
The field-induced pair creation process can be divided 
into three regimes by means of the so-called laser inten- 
sity parameter 



m e c 



(1) 



a being the amplitude of the laser potential. One dis- 
tinguishes the tunnel regime, where £ 3> 1 (and aw/c <C 
E cr ) and the pair creation is similar to a tunnelling pro- 
cess from the negative- to the positive-energy contin- 
uum, the above-threshold regime similar to the above- 
threshold ionization of atoms or molecules in strong laser 
fields where £ *~ 1, and the multiphoton regime with 
( < 1. The typical photon orders yielding the main 
contribution to the total pair production rate in the tun- 
nel and above-threshold regime are much larger than 
the minimum order ro which is required by energy- 
momentum conservation (see Eq. @ below). In the 
multiphoton regime, the main contribution to the pair 
production process comes from the minimum number ro 
allowing the creation to take place. In this regime, the 
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FIG. 1: Feynman diagram of the considered process in the 
lowest order in olj. The dashed line represents the virtual 
photon propagated between the projectile and the electron- 
positron vertices, and the double lines stand for the exact 
lepton wave functions in the laser field (Volkov states). 




FIG. 2: Decomposition of the Feynman graph in Fig. [TJinto 
a perturbation series. The left-hand side represents the total 
Feynman diagram, where the leptons are represented by exact 
lepton wave functions in the laser field, i.e. the Volkov states 
(double lines) (Eq. (O). On the right-hand side, we show a 
few examples of which underlying graphs constitute the left- 
hand side Feynman diagram. The wiggled lines represent real 
laser photons absorbed in the pair creation process. The first 
two rows on the right-hand side show the four leading graphs 
for one-photon processes, whereas the third row shows two 
exemplary two-photon processes. Higher orders involve both 
absorption and emission of laser photons, where the net ab- 
sorbed energy must surmount the energy threshold (Eq. Q). 

total pair creation rate R scales with R oc £ 2r ° . 

In the present paper, we study Bethe-Heitler e + e~ cre- 
ation in the collision of a relativistic muon with a high- 
frequency circularly polarized laser beam within the mul- 
tiphoton regime. In comparison with already existing 



calculations for projectile nuclei, which are treated as in- 
finitely heavy particles providing an external Coulomb 
field, the present approach has two advantages: firstly, 
the fact that the considered projectile particles are lep- 
tons and as such fundamental particles, allows us to treat 
the process exactly in terms of (laser-dressed) quantum 
electrodynamics up to the leading order in the coupling 
constant a/ (Fig. [I}. Secondly, the treatment of the 
projectile nuclei as external fields does not account for 
their finite mass and thus neglects any recoil effects on 
the projectile. The approach pursued in this paper al- 
lows to study the said recoil effects and their dependence 
on the projectile mass and the photon order. We show 
in particular that the muon recoil becomes the more pro- 
nounced the more laser photons (of equal total energy) 
participate in the pair production. 

The laser-dressed Feynman graph in Fig. [T] can be ex- 
panded into different orders of absorbed laser photons 
(Fig. [2|) . The pair creation process consists of all possible 
r-photon processes, which themselves are composed of all 
possible processes involving combinations of absorbed or 
emitted laser photons at each of the two vertices leading 
to r net absorbed photons. Note that, since we study 
low laser intensities, any of these diagrams could be cal- 
culated by means of ordinary quantum electrodynamics 
(QED) within r-th order perturbation theory in the pho- 
ton field. However, the framework of laser-dressed QED 
employing relativistic Volkov states is more convenient 
for our purposes [HI, [H, [H, |32j . 

In the rest frame of the incoming projectile, the photon 
energy is higher than in the laboratory system by the 
Doppler factor ^(1 + /3)/(l — (3), where (3 is the ratio 
of the muon's lab-frame velocity to the light velocity. 
Due to the finite mass of the muon projectile, not only 
the energy 2m e c 2 for the creation of an electron and a 
positron must be surmounted, but also the recoil energy 
of the scattered projectile must be provided. By calcu- 
lating the Mandelstam variable for the invariant mass of 
the process, one finds the threshold relation in the rest 
frame of the incoming muon for r absorbed laser photons 
of energy hu>, 

r ■ Huj > 2m e c 2 (l + ^) , (2) 

M being the mass of the projectile. The muon rest mass 
is M = 105.6MeV/c 2 , i. e. about 200 times the rest 
mass of an electron. Since the probability for Compton 
scattering is inversely proportional to the scattering 
particle's mass, the Compton channel occuring in the 
SLAC experiment will be highly suppressed for muon 
projectiles. 

We note that the process /i + ru — > /i + e + e~ is formally 
related to some other processes via crossing symmetry: 
for instance, electron-muon scattering in strong laser 
fields e~ /z _ — ► e _ /i~ (33j . or muon pair creation from 
positronium in laser fields e + e~ — > [H| (see also 

[IE US H3)- Moreover, there are certain similarities 
to laser-assisted M0ller scattering e~e~ — > e~e~ 13811 
and laser-assisted Bhabha scattering e + e~ — > e + e~ [39j . 
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Also, a laser-assisted Bethe-Heitler process, where a 
high-frequency 7-photon and a nucleus collide within a 
background laser field, Z + 7 — > Z + e + e~ , has been dis- 
cussed [lei]. Besides the creation of matter-antimatter 
particle pairs, there are also other interesting QED 
effects that occur in combined laser and Coulomb fields, 
such as photon fusion [4l[ or Delbriick scattering 
Efficient pair production by the linear Bethe-Heitler ef- 
fect through bremsstrahlung photons in a laser-induced 
plasma has been observed recently [43| . 

The paper is organized as follows. In Section [H] we de- 
rive an expression for the fully differential pair produc- 
tion rate. In Section IIII| the results of the numerical 
calculation of pair creation processes involving the net 
absorption of r = 1,2 or 3 photons are given, including 
total and differential pair creation rates. Section HVl gives 
a summary of the paper. In the Appendix fAl we summa- 
rize the calculation for proton projectiles and compare it 
to our calculation for muon projectiles. 
We employ a natural units system in which h = c = 1 and 
e = ^fuj {otf — 1/137 being the fine structure constant). 

We make use of the metric tensor g^ v =diag(H ), so 

that the scalar product of two four-vectors p^ — (p°,p) 
and g M = (q°,q) reads (pq) — p°q° — pq. Furthermore, 
we employ Feynman slash notation for four-products of 
four- vectors with the Dirac matrices 7**, j/= (jp). 



II. THEORY 

The Feynman graph corresponding to the process under 
consideration is shown in Fig. [TJ In the collision with the 
laser beam, the muon experiences a certain recoil and 
emits a virtual photon, which decays into an electron- 
positron pair. The amplitude of the process then reads 



d xd yif} p _ s _{x)^tp p+SJt {x) 

xV^(x~y)^p rsr {yh v ^ PS (y). (3) 



Here, x and y denote the space-time coordinates of the 
produced pair and the scattering particle, respectively. 
T^iw (x — y) is the propagator for the virtual photon prop- 
agating between the two vertices. The muon vertex is de- 
scribed by ^p'S'(y)l" , ^Ps(y), and the electron-positron 
vertex by ip ps _ (x)7 M Vv+s+ (x) with the Dirac matrices 
7 M . -0 and \t are the Volkov states [3|4^| of the respec- 
tive particles which solve the Dirac equation for spin- 1 /^ 
particles in an electromagnetic field, 



m)ip = , 



(4) 



with the laser potential = (0, A) in the radiation 
gauge. We consider a circularly polarized laser field, so 
that the four potential reads A fl (-q) — ai cos(^)+a2 sin^) 
with the laser phase 77 := kx, where k is the wave vector of 
the laser field. The four-vectors ai/ 2 read ai — (0, a, 0, 0) 
and ci2 = (0,0, a, 0), a being the amplitude of the laser 



potential. The state vector for an electron with the ki- 
netic momentum p_ and the spin projection s_ thus can 
be written as 



A>_ s _(z) = 



with the action 



Vq°_ 



1 - 



2{k P -) 



(5) 



S = -(q-x)-\ — r- sm(r)) j- r- cos(?7) . (6) 



{hp-) 



{hp-) 



Here, m = m e is the electron mass, V a normalization 
volume, and 1lp_ s _ IS £l free Dirac spinor (46[. g_ is the 
effective momentum of the electron in the laser field [43 , 



q- 



V- 



2{kp-) 



(7) 



The corresponding effective mass is = m(l + £ 2 ) with 
the laser intensity parameter £ of Eq. (fT]). The Volkov 
states for the positron can be obtained from ([5]) and ([6]) 
by replacing p^ by —p^ and u p s _ by a corresponding 
antiparticle spinor v p+s+ . Replacement of the coordinate 
7] = kx by k = ky, the mass m by the projectile mass 
M, the kinetical momentum p_ by P or P' for the in- 
coming or scattered muon, respectively, and the effective 
momentum g_ by Q or Q' yields the Volkov states for 
the initial and scattered muon. The effective muon mass 
is obtained by M* = M(l + S 2 ), and the corresponding 
intensity parameter reads S = ea/M — £m/M. 
We concentrate our study on the multiphoton regime 
£ <C 1, where the intensity of the laser field is com- 
paratively low and the photon energy is high. Because 
of the low intensity, the electric field is much smaller 
than the critical value E cr = 1.3 x 10 16 V/cm 0. There- 
fore, we neglect vacuum polarization effects which be- 
come important only at near-critical laser intensities 
I > 10 29 W/cm 2 [43, |H, SI] and employ a free pho- 
ton propagator for the description of the virtual photon 
which is propagated between the two vertices: 



y) 



d 4 q ^ei^-v) 
(2?r)* g 2 
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(8) 



with the integration variable q describing the momentum 
of the virtual photon. 

It is important to note that the amplitude in Eq. ([3|) 
fully accounts for the interaction of the leptons with the 
laser field by employing the exact solutions to the Dirac 
equation, i. e. the Volkov states; the interaction between 
the leptons and the QED vacuum however is taken into 
account only to lowest order in qj. A similar approach 
can be found in 0, M, H ^ . 

The space-time integrals in ([3J can be performed by ex- 
panding the individual vertex expressions into Fourier 
series. Using the generating function of the regular cylin- 
drical Bessel functions [50( , one finds 



4 



a f 2(2nfMm 
~i V 2 



d 4 q 



q 2 sJq° + q°-Q' Q n,N 



J2 M^ie+e- \n)M„(n, n'\N)6(q + q+ + q-- nk)5(q + Q + Nk-Q'). (9) 



The integer numbers n and N correspond to the numbers 
of absorbed laser photons at the electron-positron and 
projectile vertex, respectively. Note that n and N may 
become negative, amounting to photon emission at the 
respective vertex, whereas the total number of absorbed 
photons, r = n + N , may not. The electronic spinor- 
matrix product A4 ll (e + e~ \n) can be written as 



, e \n := u 



P—S- 



X 



2(kp+)(kp_) 



+ ^W7)'' w ' ~ Jik)^ ' D ")K" + • (10) 

with the coefficients 
B n = J„(o)e' m "° , 

C n = i(j n+ i(o)e i ( n+1 i*» + Jn-iCaJe*^*) , 

D n = ^(j n+ i(a)e^ n+1 ^ - J^x^e^"- 1 '" ) . (ff) 

The functions Ji(a) are the regular cylindrical Bessel 
functions of integer order. Their argument is a — 
\A*i + a 2 an< i the angle 770 is given by cos 770 = OL\j a, 
sin 770 = 0:2/0: with 



(hp-) (kp+) 



J = 1,2. 



(12) 



Note that the Bessel functions account for the fact that 
in the multiphoton regime, the main contribution to the 
total pair production rate comes from the minimum num- 
ber of absorbed photons r*o, for which the threshold re- 
lation ^ is fulfilled: for small arguments a < 1, the 



Bessel functions scale with J„(o) 



a 



n\ 



[50|. Since the 



laser field parameter £ is supposed to be small and it can 
be factored out in the argument of the Bessel functions, 
higher photon orders yield partial rates which are several 
orders of magnitude smaller than the one for ro. 
The corresponding muonic spinor-matrix product 
M„((J,,n'\N) is 



M v 0*,MW :=EW (7" + 



2(kP')(kP) 

+ ^j^^ + j^^' Hn )) Ups - (13) 



r 



The different signs in Eq. (JTUJ) and (fT5|) occur because the 
former spinor-matrix product describes a particle and 
an antiparticle, whereas in the latter only one type of 
particle appears. The argument of the Bessel functions 
Jn(/3) for this vertex reads 



(14) 



with 



{kP>) 



(kP) 



(15) 



Making use of the angle k defined by cosk = fli/fl, 
sinK = the coefficients Fjf, Gn, and Hm are de- 

fined similarly to Eq. (fTTj) involving Jn((3). 
We can separate the amplitude into a sum of partial am- 
plitudes for one particular photon order r = n + TV, 



&<r) =M'Y^M' l (e + e-\n) ■ M M (/J, fi'\r-n) 

n 

5{q+ +q _ + Q'-Q- rk) 
(q + q + + q- - nk) 2 

with the factor 



jV y = £/_ 2(27r) 5 Afm 



(16) 



V 2 Jq + q°_Q' 0n0 



(17) 



Then, the total transition amplitude and its square read 

y=^y {r \ \y\ 2 =^2\y {r) \ 2 , (is) 

r>ro r>ro 



where ro is the minimal number for which the threshold 
relation {2J is fulfilled. Because of the delta function in 
Eq. p^|) , there is no double sum over r, r' so that the 
summation in (|18p can be pulled out. The square of the 
partial amplitude y^ is 



5 



|^(r)|2 = |AA'| 2 ^ 



■M^(e+e-|n).M M (//, //|r - n)A)t( ft //jr - nQA^e+e- K)?V% + + g _ + g/ _ q _ rfc ) 

(27r) 4 (g+ + g_ - nfc) 2 ■ (g+ + g_ - n'fc) 2 



(19) 



The factor TV/ (2tt 4 ) stems from the square of the delta 
function [4r| . leaving a time factor T and a volume fac- 
tor V which is the same as in the normalization of the 
Volkov states flS)). 

The total pair production rate is obtained from (|19p by 
averaging over the possible initial spin states, summing 
over the final spin states, integrating over the final mo- 
menta, and dividing the result by a unit time T: 



R = 



Vd 3 Q' f Vd 3 q+ f Vd 3 q- 1 

J l^y J ~W?2 



^E E 

S S f ,8+,S- 



T 



(20) 

With the completeness relation for the components of the 
Dirac spinors u p _ s _ and v p+s+ [46T ]. we can perform the 
spin summation by introducing the matrices 



for the projectile vertex and 



r L n 



(- 



e 2 a 2 k^ 



kp. 



(22) 



A™ — 



7/i 



e 2 a 2 k tJl k' 



F r 



i 



2 VfcP 



1 

fcp 



(21) 



for the vertex of the produced pair. We find 



M^e+e-^M^, p!\r - n)Ml{p, n'\r - n')M^{e + e -\n') 

S,S',s+ ,s- 



Tr 



( r M - 



m ™ i- + m " 



i -) • i: (..a:: a:: — — ) = : r: n 



Ira r n 2m 



2M r v 2M 



(23) 



This trace product can be evaluated by the standard pro- 
cedure. 

Since the kinematics of this process are rather compli- 
cated in the rest frame of the incoming muon, we now 
transform into the center-of-mass (cm.) system in order 
to perform the integrations in Eq. (|20|) . We can separate 
the total rate into partial rates R = ^2 r R r correspond- 
ing to the contributions of particular photon orders r like 
we did for the transition amplitude. For each photon or- 
der r there is a corresponding cm. system. Let u> be the 
laser frequency in the rest frame of the incoming projec- 
tile muon and k — u>(l, 0,0,1) the corresponding wave 
vector. In the cm. system, the incoming muon moves 



along the z-direction and 

[QlmA ( 

Q cm ,r = Qln,r = ~T ■ Km,r = -r \ ) , (24) 
\Qcm,rJ \Wcm,r , 



where the index (cm, r) indicates the cm. system cor- 
responding to the photon order r. The laser frequency 
transforms as 



rUlnm.r = TUJ A 



with the velocity 



! 1-(3 C 



1 ~t~ ftcm.r 



M* + 2rui 



0C 



IQc 



ruj 



>> crn.r 



(25) 



(26) 
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electron, 



FIG. 3: (Color online.) Visualization of the coordinate trans- 
formation which Eq. (|27p is based upon. 



According to this, the center of mass as viewed from 
the projectile rest frame moves with the Lorentz fac- 
tor 7 r = (1 — /3^ m r ) -1 / 2 . We can perfom one of the 
momentum integrations in Eq. (|20[) by making use 
of the three-dimensional 5-function. This gives e. g. 
(q + )em,r = ™Qcm,r ~ (q_)cm,r< For the sake of brevity, 
we will omit the index (cm, r) where there is no ambigu- 
ity. From now on, all quantities refer to the cm. system. 
The remaining integrations are performed in spherical 
coordinates. Because of the spherical symmetry we may 
choose the x-axis of the coordinate system freely. Thus, 
we may set the azimuth angle of the scattered muon 
4>' = so that the integration yields the constant fac- 
tor 2ir. 

Regarding the remaining integrations, we perform a co- 
ordinate transformation following Mork [5l| . We regard 
the solid angle element of the created electron within a 
coordinate system where the z'-direction is given by the 
momentum vector of the scattered muon (see Fig. [3]). 
For this transformation, the Jacobi determinant is unity, 
J = 1, so that the partial rate in the cm. system reads 



R' 



a 2 f m 2 M 2 



/7 |q-HQ'l V- T „„' 

J v Q° + Q° ' V r 



2tt 2 



5(E-q°, 



q°- 



') 



(<7 + +q- — nk) 2 (q + + q_ — n'k) 2 
■ d(cos#'_)d(p'_d(cose')dq°_dQ' 



-=-Q'-q_ 

(27) 

with the abbreviation E := Q° + rk°. The new polar 
angle of the electron, , is the angle between the mo- 
mentum vectors of the electron and scattered muon and 
is found in the energy-conserving (5-function (because 
q + = Q' q ). The argument of the 5-function has 
one root at 



cosC 



(£-Q'°) 2 -|QT-2g° (E-Q'°) 
2|Q'||q_| 



(28) 



and from the requirement that |cosi?'_!?| < 1 follow the 
integration limits for the energy q°_ of the created 



Q'°±\Q'U l- 



2m? 



E{Q0-Q'°) 



(29) 



Since these limits must be real numbers, one finds the 
upper limit for the integration over the energy Q'° of the 
scattered muon, 



M 2 _ 

mi 



(30) 



and the lower limit is M*. 

Thus we can write the partial rate in the cm. system 

f+l r 2ir 



R, 



dQ 



HI 



dq°_ 



d cos 9' 



d^'_Af 



E 

n,n 



rpr< 



\ (q + + q- — nk) 2 {q + + q~ — n'k) 2 



with the factor 



a^m 2 M 2 |q_||Q'| 



COS —COS 1?^ 

q + =-Q -q_ 

(31) 
(32) 



2^ 2 ■ q° + Q° ' 

Eq. (f3"Tj) gives an expression for the partial rate cor- 
responding to one particular photon order r in the re- 
spective cm. system. In order to obtain the total pair 
creation rate, we have to sum over all relevant orders in 
the rest frame of the initial projectile: 



r=r r=r 



(33) 



where the Lorentz factor "f r accounts for the time dilation 
between the frames. In the multiphoton regime consid- 
ered in the following, all the summands with r > tq yield 
negligibly small contributions, so that only the leading 
term with ro absorbed photons must be taken into ac- 
count. The rate (f3"3"|) may be transformed into the labo- 
ratory system according to 

i? Lab = — , (34) 

7Lab 

where 7 Lab is the Lorentz factor of the incoming muon in 
the laboratory system. Note that total probabilities of 
multiphoton processes are usually expressed directly as 
rates rather than cross sections, because the cross section 
ay = Rr/j of a nonlinear process (r > 1) still depends on 
the incoming photon flux j cx £ 2 , which is an undesired 
feature for a cross section in the usual sense. 



III. RESULTS 
A. Total pair creation rates 

In this section, we present the results of our numerical 
calculations of the total pair production rates. They re- 
fer to the rates of produced electron-positron pairs per 
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projectile muon in the restframe of the incoming muon 
for a laser beam of infinite length. In order to obtain the 
total pair yields in the laboratory system, which could 
be observed in a corresponding experiment, one has to 
transform into the laboratory frame (see Eq. (j34|0 and 
multiply the outcoming rate by the interaction time. The 
latter is given by half the laser pulse duration r, provided 
that this is considerably small compared to the muon 
lifetime, which usually is the case (typically, r ~ fs — ns, 
while the muon lifetime in the lab frame is r„ = 7 Lab -2/is). 
Finally, one has to multiply by the number of muons in 
the projectile beam. 

1. Linear process 

For linear processes, the corresponding expressions for 
the total rates can be found using the framework of stan- 
dard QED jH |H [E|. This allows for a testing of our 
calculation by comparison. 

We first consider an e + e~ pair creation process by a sin- 
gle photon of the energy to — 1.8MeV impinging on a 
muon initially at rest. We choose the laser intensity pa- 
rameter to be £ = 7.5 x 1CP 4 , which corresponds to XFEL 
photons with a lab frame energy of w Lab = 9keV and an 
intensity of I Lab = 8 x 10 19 W/cm 2 , which are envisaged 
to become available by the X-ray lasers planned at SLAC 
(Stanford, USA) and DESY (Hamburg, Germany) 
In this constellation, a muon Lorentz factor of 7 Lab = 100 
would be necessary. Alternatively, if even higher muon 
energies were to become accessible [H, [H, [13], the as- 
sumed set of parameters could be implemented by an 
XUV laser with a photon energy of w Lab = 90eV and 
an intensity of 7 Lab = 8 x 10 15 W/cm 2 (58| with a muon 
Lorentz factor of 7 Lab = 10 4 . The numerical calculation 
of the total rate in the rest frame of the incoming muon 
(f2T))) yields for these parameters 

Rr=i(u = 1.8MeV) = 1.28 x 10V 1 . (35) 

This rate results from a single term in the double sum 
of Eq. (|3~Tj) with r = n = n' = 1, i. e. one laser photon is 
absorbed at the electron-positron vertex, and no photon 
is absorbed (or emitted) at the projectile vertex. This 
term comprises the first two Feynman diagrams on the 
right-hand side in Fig. [5] All other diagrams yield 
negligibly small contributions. 

The rate is proportional to £ 2 . In mathematical terms 
this scaling arises from the Bessel function J n (a) in 
Eq. (JTUJ). Since a <C 1 and a oc £, the amplitude for 
one-photon absorption at the e + e~ vertex is propor- 
tional to £, leading to R r= \ cx £ 2 . The properties of the 
Bessel functions also explain why terms involving photon 
exchange at the muon vertex are strongly suppressed 
and may be ignored. Those terms are proportional 
to J N {j3) oc ^ |JV| with (3 oc S = £m/M, leading to 
contributions to the process rate which are smaller than 
the leading term (N — 0) by at least three orders of 
magnitude. Higher photon orders yield even smaller 



contributions. 

Since the recoil energy is small compared with Q° sa M, 
the main contribution comes from the [i = v = term 
in the sum over the spinor-matrix product (|10[) (see 
also App. [A}. We found this term to yield 97.5% of the 
given result (pj)) . 

We expect the numerical value of the total pair pro- 
duction rate for muon projectiles to be approximately 
the same as for e. g. proton projectiles because both 
particles are very heavy as compared to the electron 
mass. Since protons are not fundamental particles, our 
calculation does not actually apply to them. But since 
they are spin-Y2 particles like leptons, we may approxi- 
mately treat them as effective Dirac particles and check 
the result of our numerical program code for projectiles 
having the mass of a proton. This yields a total pair 
production rate of =m "(1.8MeV) = 1.31 x 10 9 s~\ 
which is in agreement with the result obtained from the 
calculation treating protons as external Coulomb fields 
[28| . The difference between the total rate for muon and 
proton projectiles stems from the recoil effects due to 
the lower mass of the former particle. 

If the absorbed energy is very close to the threshold en- 
ergy ©, i. e. the photon energy in an r-photon process 
is close to 



one finds a power law for the dependence of the total pair 
creation rate on the amount of energy absorbed addition- 
ally to u) r mjJL : 

R r oc (lj — w r . min ) Xr . (37) 

For one-photon e + e~ pair production processes, our nu- 
merical results for muon projectiles can be fitted to the 
functional relation of Eq. (|3"T|) and yield the exponent 

Xi - 2.97, (38) 

with an accuracy of 0.1% within the range between ix>i ;min 
and 1.05 x L0i imiu . In the case of very heavy projectile 
particles, the exponent is xf 1 ^ 00 = 3 [Hj]. An investi- 
gation of the scaling behavior for proton projectiles with 
our program code yielded agreement with this value. For 
electron projectiles, the exponent is xT" = 2 |53|- We do 
not consider electron projectiles here, since our calcula- 
tion does not take into account the quantum mechanical 
exchange term of the produced and scattered electron. 
But we may calculate the rates for hypothetical pro- 
jectile particles carrying e. g. twice the electron mass. 
While the minimum absorbed energy wi jmin is 4m e for 
electron projectiles, and approximately 2m e for protons 
and also for muons, for hypothetical projectile particles 
with M — 2m e , the threshold energy lies just in the mid- 
dle at 3m e . For such particles, it is not possible to neglect 
summands with N, N' ^ in ([31]) . Instead, we employed 
the two leading orders, n = n' = 1 and n = n' = 0. We 
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find the exponent to be " 2, " e = 2.02 with an accuracy 
of 0.2% between w l min and 1.06 x uj l nlin . This exponent 
is close to the one obtained for electron projectiles. 

2. Nonlinear processes 

The same amount of energy as in the one-photon process 
considered in the foregoing paragraph is absorbed in a 
two-photon process with half the photon energy, i. e. 
lu = 900keV in the incoming muon's rest frame. Using 
the same laser intensity parameter as in the linear case, 
we obtain the total pair production rate for the two- 
photon process with muon projectiles, 

i?r= 2 (w = 900keV) = 187.5s" 1 . (39) 

This is considerably less than obtained for the absorption 
of one single laser photon with twice the energy. Since 
the pair production rate scales with £ 2r , we expect the 
rate for the two-photon process to be roughly 7 orders 
of magnitude smaller than for the one-photon process, 
which agrees with Eqs. (|3"5|) and (f3"9"|) . Also here, the 
main contribution to the pair production rate comes 
from the fi = v = summand in the trace product ([23]) 
according to (|A12p . In this case the said summand yields 
97.8% of the total rate (f3"9")) . The calculation neglecting 
the finite proton mass yields for this parameter constel- 
lation the value of R^°°(u = 900keV) = 190.0s" 1 

As in the foregoing paragraph, it is possible to find a 
scaling behavior of the total pair production rate for total 
absorbed energies close to the threshold energy, see Eq. 
(|37p . For two-photon processes we find the exponent 

X2 = 3.92 (40) 

for muon projectiles with a relative error of 0.1% between 
<^2,min and 1.05 x W2,min- Considering proton projectiles, 
we find the exponent to be 3.96 with an error of 0.23%, 
which agrees with the result in [28| , and for hypothetical 
particles with twice the electron mass, the exponent is 
3.02 with the accuracy of 0.3%. 

We observe that in each case, the exponents exceed the 
ones found in the foregoing paragraph by about 1. An 
investigation of higher photon orders up to r = 5 shows 
that the increase of the exponent is proportional to the 
increase of the photon order (see Fig. H}. We find 

Xr»1.95 + r (41) 

for muon projectiles. 

It should be noted that regarding this scaling behavior, 
the circular polarization of the laser is crucial; for lin- 
ear polarization, the scaling would be different in nonlin- 
ear processes. There one finds x^ 7 ^ 00 = 2, for instance 
[28| . Polarization dependence is a characteristic feature 
of nonlinear processes, in general. 




FIG. 4: Pair creation rates R r for muon projectiles as a func- 
tion of the energy absorbed additionally to u Timin of Eq. 1)360 
for photon orders r = 1,2,3. The exponent increases as the 
number of absorbed photons increases. 

As in the one-photon process considered above, only 
N = N' = summands significantly contribute to the 
total pair creation rate, i. e. the absorption or emission 
of laser photons at the muon vertex are suppressed. We 
note however, that so-called resonances can occur in the 
terms involving photon exchange both at the electron- 
positron vertex and at the muon vertex (n ^ 0,N =/= 0). 
Then, for certain final state momenta, the virtual photon 
may reach the mass shell, q 2 = (q+ + q- — nk) 2 = 0, and 
become a real photon. Physically, this resonant behavior 
means that the original second-order Feynman graph in 
Fig. [1] decomposes into two first-order processes: first 
the emission of a real photon by a muon and second the 
production of an e + e~-pair by a real (non-laser) photon 
(Breit- Wheeler process). Both processes are kinemati- 
cally allowed in the presence of an external laser field; the 
second process was responsible for the pair production 
observed in the SLAC experiment [U HI] • In our case 
however, the real-photon channel is strongly suppressed 
by the large muon mass: while the process involving the 
propagation of a virtual photon with no absorption at 
the muon vertex yields a production rate proportional to 
£ 4 for a two-photon process, the corresponding rate for 
pair production involving the absorption of one photon 
at the muon vertex is proportional to £ 2 E! 2 <C £ 4 . Res- 
onances in multiphoton processes have mainly been dis- 
cussed for the examples of laser-assisted e~e~ (M0ller) 
scattering [38[, e + e~ (Bhabha) scattering [3i|, and elec- 
tron bremsstrahlung [49(. 



B. Recoil effects 

In this section, we discuss the recoil distributions in 
the cm. system for the already considered one- and 
two-photon processes, and also for a three-photon 
process, for different projectile particles. 
Figure [5] shows the cm. differential rate for muon 
and proton projectiles plotted against the energy loss 
relative to the incoming projectile's kinetic energy. 
Note that for the protons, we divided the differential 
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FIG. 5: (Color online.) Differential pair creation rate plotted 
against the energy loss of the projectile relative to its initial 
kinetic energy for the one-photon process at u — 1.8MeV. 
The red solid line shows the recoil distribution for muons, the 
green dotted curve shows the one for proton projectiles. For 
protons, the differential rate has been divided by a factor 10. 



rate by a factor 10 in order to fit into the graph. 
In the cm. system, the incoming muon projectile 
has a kinetic energy of Ekin = 14.83keV, and after 
the collision with the laser beam, the most probable 
kinetic energy is E' kin — 8.81keV, i. e. the muon loses 
about 40.6% of its initial energy of motion. The 
proton projectile possesses an initial kinetic energy of 
Ekin = 1.72keV in the cm. system, and E' kin = 1.03keV 
after the collision; i. e. it loses about 40.2% of its 
initial energy of motion. The absolute loss of energy 
is about Ekin — E' kin =: AQ° — 6keV for muons and 
AQ° = 0.7keV for protons, which is in agreement with 
the general rule that the recoil momentum A|Q| should 
be of the order of 2m. Since the heavy projectiles move 
nonrelativistically, this should lead to a recoil energy 
AQ° ~ A\Q\ 2 /2M, and therefore the absolute energy 
loss for the muon should be one order of magnitude 
larger than for the proton. As discussed in Section 
IIII Al and in Appendix [A] the total pair creation rates 
for muon and proton projectiles yield almost the same 
value. As shown in Fig. [5j the shapes of the recoil 
distribution for the two projectile species are similar. 
The distribution for the proton is about a factor 10 
higher than the one for muons, but it is also a factor 10 
narrower if absolute energy loss is considered, so that 
the integration leads to about the same numerical value 
for both projectile types. 

If the same amount of energy is absorbed in a two- 
photon process, as considered above, the shape of the 
recoil distribution is quite similar to the one for the 
one-photon process in Fig. [5J but slightly broader and 
shifted to higher relative energy loss (Fig. 0. The 
initial kinetic energy of the projectile particles in the 
cm. system is the same as in the one-photon process, 
since the total absorbed energy is the same in both 
cases. Here, for the muon the most probable kinetic 
energy after the impact of the laser beam is 7.93keV 
which is tantamount to a loss of 46.5% of the initial 
kinetic energy. The proton projectile loses 46.1% of its 
initial kinetic energy and ends up with a kinetic energy 
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FIG. 6: (Color online.) Same as Fig. but for a two-photon 
process with uo = 900keV. 



of 0.93keV. 

In both the linear and the nonlinear process, both 
projectile species lose about the same ratio of their 
initial energy of motion in the collision with the laser 
photons. But although the total amount of absorbed 
energy, as well as the initial energy of the projectile 
in the cm. system and the upper limit (|30[) for its 
final energy are the same in both considered processes, 
the lost energy is significantly larger in the nonlinear, 
two-photon process. 

This trend is confirmed by the investigation of another 
kincmatically identical process yielding the same value of 
the totally absorbed energy from the laser field: if three 
photons of energy 600keV are absorbed, the scattered 
muon's most probable kinetic energy is 7.49keV, which 
corresponds to a loss of 49.5% of the initial kinetic 
energy. The most probable kinetic energy of a scattered 
proton projectile lies at 0.88keV, which means a loss of 
49.1% of the initial kinetic energy. It should be noted 
that the considered three-photon process is unlikely to 
be experimentally observed, because for the considered 
photon energy, also a two-photon process is possible. 
Because R r oc £ 2r , the two-photon process is 7 orders 
of magnitude more probable than the three-photon 
process. 

We point out that the increasing of the energy loss 
of the projectile is in agreement with the fact that 
the emission angles of the produced pair rise with the 
photon order [26j. The angular distribution of the 
electron and positron is mainly governed by the Bessel 
function J n (a), where a oc sin?9 9± (see Eq. (|12[) ). Due 
to the small-argument behavior of J n (a), large values 
of n favor large emission angles i? g± . Since the typical 
energy of the electron and positron is not much affected 
by the photon order (see Fig. [8] below), this implies 
that the longitudinal component of the scattered muon 
momentum decreases when n (i. e. r) increases (see also 
Fig. [TU] below). 

In order to demonstrate that our numerical calculation 
yields larger recoil effects for lighter particles, we consider 
once more the case of a hypothetical particle possessing 
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FIG. 7: (Color online.) Recoil distribution for hypotheti- 
cal projectile particle with mass M = 2m for a one-photon 
process at to — 1.8MeV. The two most important orders, 
n = n' = 1 (red solid line) and n — n' — (green dotted 
line) are depicted. 



twice the electron mass. Figure [7] shows the recoil distri- 
bution for the main contributions to the total rate stem- 
ming from n — n' = 1 and n — n' = for the case of a 
one-photon process at the photon energy ui = 1.8MeV in 
the restframe of the incoming projectile. Initially, the hy- 
pothetical particle possesses a kinetic energy of 305keV. 
The total recoil distribution consists of the sum of the 
two depicted orders and has its maximum at 263keV, 
i. e. the particle loses 86.4% of its initial kinetic energy. 
As expected, the lighter particle loses a higher ratio of its 
initial energy in the collision with the laser beam. Also 
here, our calculation of the corresponding two-photon 
process shows that the energy loss is larger in the non- 
linear than in the linear process. 



C. Differential rates 

In this paragraph, we discuss some further differential 
pair production rates for muon projectiles. Each of the 
following figures is to be understood in the cm. sys- 
tem for the respective r-photon process. As in the fore- 
going paragraph, we consider a one-photon process at 
u; = 1.8MeV, a two-photon process at oj = 900keV, and 
a three-photon process at uj — 600keV. 

Fig. [5] shows the dependency of the (normalized) differ- 
ential pair creation rate on the energy of the produced 
electron. The most probable electron energy, at which 
the maximum occurs, is about 890keV for all considered 
photon orders. It slightly decreases with increasing pho- 
ton order, but the difference between the most probable 
electron energy for the one- and the three-photon pro- 
cess is only 0.7keV. 

The differential rates in Fig. [5] have been divided by 
their respective maximum values. This allows for a com- 
parison of the distributions for the different photon or- 
ders. The curvature of the distribution increases with 
the number of absorbed laser photons: for the consid- 
ered one-photon process at w = 1.8MeV, the relative 



FIG. 8: (Color online.) (Normalized) double differential pair 
creation rate in arbitrary units as function of the energy of 
the produced electron for muon projectiles. The thick red line 
shows the distribution for the one-photon process at 1.8MeV, 
the green dotted line shows it for the two-photon process at 
900keV, and the thin blue line for the three-photon process at 
600keV. The kinetic energy of the scattered muon is set to the 
most probable value, i. e. 8.8keV for the one-photon, 7.9keV 
for the two-photon, and 7.5keV for the three-photon process. 
The distributions have been divided by their maximum values 
in order to show the evolution of the shape with increasing 
photon order. 




FIG. 9: (Color online.) Dependency of the differential pair 
creation rate on the scattered muon's polar angle. The respec- 
tive photon energies and the corresponding scattered muon's 
energies are as in Fig. [5] The energy of the produced 
electron has been set to the most probable values, which is 
q°_ « 890keV in all cases. 



height of the maximum is only 0.4% of its absolute value. 
Considering the two-photon process at u) = 900keV, we 
find a similar shape, but here the relative height of the 
maximum is 15% of its value, and for the three-photon 
process, the relative height is 37%. The reason is that 
in processes of high photon order, the e + e~ pair is pref- 
erentially produced in a symmetric configuration with 
q° + « q°_ w (E - Q'°)/2 and tf+ m (see Ref. (H). 
Therefore, deviations from the maximum in the middle 
of the allowed energy range yield smaller contributions 
and the maximum becomes more pronounced. 

Fig. [9] shows the dependency of the differential pair pro- 
duction rate on the cosine of the scattered muon's polar 
angle for the considered one-, two-, and three-photon 
processes. The higher the number of absorbed laser pho- 
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FIG. 10: Visualization of the kinematics in the c.m. sys- 
tem. The r laser photons are depicted as one "big" photon. 
The vanishing total momentum in this frame leads to a cor- 
responding correlation of the particle momenta in the final 
state. 

tons, the less the muon is deflected from its initial di- 
rection of motion along the negative z-direction. (Note 
that the laser beam propagates in positive z-direction, 
see Fig. [TO]) This behavior is best understood by consid- 
ering the argument used to explain the increasing curva- 
ture of the produced electron's energy distribution with 
increasing photon order. Since the symmetry between 
the produced electron and positron emission angles be- 
comes more pronounced for higher photon orders, the 
transverse momentum of the scattered muon and thus 
its deflection angle are reduced (see Fig. fT0|) . 
Fig. [TUJcan also help us visualize the aforementioned in- 
crease of the muon's energy loss. It has been said that 
for large photon orders n = r, the Bessel functions J n {a) 
favor the emission of the e + e~ pair under large emission 
angles. The typical energies of the electron and positron 
are unaffected by the photon order (Fig. [8]). Since the 
emission angles of the electron and positron get larger 
while their energies stay the same, the longitudinal com- 
ponent of the sum of the pair particles' momenta gets 
smaller for higher photon orders. Therefore, as can be 
seen in Fig. [TUl the scattered muon's longitudinal mo- 
mentum projection must decrease with increasing photon 
order. 



IV. CONCLUSIONS 

In this paper we considered the creation of electron- 
positron pairs in the head-on collision of a relativistic 
muon and a high-frequency laser beam in the multipho- 
ton regime. 

The results for the total pair creation rates for linear 
and nonlinear pair production have been compared to 
the results known for projectiles of infinite mass. Close 
to the threshold energy, we derived a power law for 
the pair creation rates' dependence on the absorbed 
energy surmounting the threshold. We found the expo- 
nent to increase linearly with the number of absorbed 
photons. The corresponding recoil distributions have 
been compared to the ones of protons as effective Dirac 
particles and to a hypothetical lighter particle of twice 



the electron mass. We have shown that in the c.m. 
system, the projectile's loss of kinetic energy in the 
inelastic collision with the laser beam relative to its 
initial kinetic energy is approximately the same for 
muon and proton projectiles. For both species, we found 
characteristic multiphoton signatures: the projectile 
energy loss increases with the number r of totally 
absorbed laser photons (with ruo = const), while the 
deflection angle of the scattered projectile decreases. 

Our results could be tested experimentally by combin- 
ing a coherent high-energy photon source with a beam 
of relativistic muons. Today, muon beams with Lorentz 
factors of about 7 Lab = 2.8 are produced e. g. at PSI 
in Switzerland or TRIUMF in Canada [59] . Usually, 
they are then slowed down in order to produce e. g. 
muonic atoms. In the present paper, we consider the low- 
intensity high-photon-energy regime of the pair produc- 
tion process. Therefore, we need high photon energies in 
the rest frame of the incoming projectile muon. Assum- 
ing a lab-frame photon energy of 9keV, as is planned to 
be reached at the upcoming XFEL facilities [5J| , a muon 
Lorentz factor of 7 Lab = 50 would be needed to obtain 
a rest frame photon energy of 900keV, the photon en- 
ergy we assumed for the considered two-photon process. 
An according muon acceleration is necessary for experi- 
mental observation of the here discussed nonlinear pair 
production processes. Corresponding efforts are being 
undertaken to develop muon accelerators for e. g. neu- 
trino factories or muon colliders reac hing muon energies 
up to the TeV range ( 7Lab ~ 10 4 ) [H,!!!^. 

APPENDIX A: COMPARISON BETWEEN 
MUON AND PROTON PROJECTILES IN THE 
MULTIPHOTON REGIME 

In this appendix we want to compare the present treat- 
ment of the nonlinear Bethe-Heitler process for muon 
projectiles to the external Coulomb field approach [T|| 
M, EH, H3, [H, [H, [H, [H, H3, EH, HI HI for proton pro- 
jectiles. To this end, we first summarize the main steps 
and results of the latter approach. 

The heavy projectile is taken into account only as the 
Coulomb field Ac of a nucleus, so that a proton in its 
rest frame can be described by the four potential 

A£ = (^,0,0,0). (Al) 

The approximation of an infinitely heavy projectile leads 
to an amplitude which describes the transition of an elec- 
tron from a negative- to a positive-energy state as in- 
duced by the Coulomb field of the proton: 

^p+p- = ie J tpp-^s-Ac^p+^+d^x , (A2) 

with the Volkov states of Eq. (T5J) . Similarly to the expan- 
sion for the electronic vertex in the calculation of Sec. [Ill 



12 



this expression is expanded in a Fourier series with coef- 
ficients depending on the regular Bessel functions J r (a). 
One finds for the partial rate corresponding to the photon 
order r 



R<r-> = 



Vd 3 q + Vd 3 q_ 1 
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with the square of the partial transition amplitude 
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(A4) 



where the four vector q r = q + + q_ — rfc represents the 
momentum transfer on the proton. Since the proton is 
assumed to be infinitely heavy, there is no energy trans- 
fer, g° 
reads 



0. The spinor-matrix product M.p/ P _ here 
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The coefficients B r , C r and D r are defined as in Eq. (fTTjl . 
Note that in this approach, all the photons are absorbed 
at the vertex of the produced pair (since the projectile 
vertex is not taken into account), r = n, independently of 
the considered laser intensity regime. Because the proton 
is treated as an external Coulomb field, only the /x = 
components of Eq. pH|) contribute to the transition. 
Thus, one finds the expression for the spin sum over the 
square of the spinor-matrix product 



+p- 



= Tr H 



2 m 



-r 



2m ' 



(A7) 



With this, we have listed the most important relations 
of the calculation for proton projectiles. We now show 



that the exact treatment for muon projectiles performed 
in Sec. HI] is structurally similar to it in the limit of small 
laser intensities. For small 5 <C £ <C 1, the effective mass 
and momentum of the muon read 



2^ 2 



M* = iliVl + S 2 w M , Q" = P" + E 



P" . 

(A8) 

Because of the large muon mass, we may assume nonrel- 
ativistic momentum changes for the scattered projectile. 
Then, in the rest frame of the incoming muon, the energy 
of the scattered muon is close to its rest mass: 



\pr\2 

P'°«M + L_L, (A9) 
with \P'\ -c M. Wc now consider the argument of the 
Bessel functions governing the muon vertex, (3 from Eq. 
(fT?]). It can be written as 



= 



co(P^-Pi) M 



< 1. 



(A10) 



Since, for small arguments, the regular Bessel functions 
scale with Jn(P) cx /S'^' [50], all orders TV ^ are neg- 
ligible for very small (3. Thus, the double sum in (|3lj) 
collapses to a sum over n — n' — r. This means that in 
the limit of nonrelativistic motion of the scattered muon 
within the rest frame of the incoming one, no emission 
or absorption of laser photons takes place at the projec- 
tile vertex. In this case, the summands in the matrix 
elements (|2"T1) estimate to be negligible except for j^Fq. 
The Bessel function Jo is about unity for small argu- 
ments, Jo(0) = 1 [50| . Then, the muon part of the spin 
sum over the spinor-matrix product in (|23p becomes 



^M / ,(/iy|0)A4(/i,//|0) 
+M 



s,s- 



M 



2M 



2M 



(All) 



so that the total sum over final spins and average over 
initial spins (f2"3")l estimates 



M^(e + e-\n)Mp(fi^'\r~n)Mi^^'\r-n')M^(e + e-\n') 

S S' : s + ,s_ 

= \ll E M%e+e-\r)M ^,^\0)MU(,,mM^e-\r)^ (A12) 

I 



which is just the expression obtained for projectiles of infinite mass, Eq. (|A7|) . 
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Another formal difference between the two approaches is 
that in the present calculation, we divide the transition 
amplitude by the four-momentum of the virtual photon 
(Tin]) , whereas in the external Coulomb field approach the 
division by the virtual photon's three-momentum is per- 
formed (|A4j) . Performing the space-time integral over 
the projectile coordinate in ([3]) under the assumption of 
no photon emission or absorption at the projectile ver- 
tex, the obtained four-dimensional <5-function (see Eq. 
©) reads 5(Q + q — Q'). With the approximation (|A9|) . 



only S(q°) remains in the 0-component of this 5-function 
(cf. Eq. IA4|) . Therefore, \q 2 \ = q 2 , and the difference 
between the two approaches concerning the division by 
the square of the virtual photon's momentum vanishes 
in the considered limit (see also [60], Sec. 4.8 for a simi- 
lar argument). In the present calculation, the remaining 
three-dimensional 5-function allows for the additional in- 
tegration over the scattered projectile's momentum to be 
carried out. 
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